Insights into molecular conduction from I-V asymmetry 
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We investigate the origin of asymmetry in the measured current-voltage (I-V) characteristics of 
molecules with no inherent spatial asymmetry. We establish that such molecules can exhibit asym- 
metric I-V characteristics due to unequal coupling with the contacts. In contrast with spatially 
asymmetric molecules, conduction takes place through essentially the same level in both bias di- 
rections. The asymmetry arises from a subtle difference in the charging effects, which can only be 
captured in a self-consistent model for molecular conduction. For HOMO-based conduction, the 
current is smaller for positive voltage on the stronger contact, while for LUMO conduction, the 
sense of asymmetry is switched. 
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PACS numbers: 85.65.+h, 73.23.-b,31.15.Ar 

Future electronic devices are quite likely to incorporate 
molecular components, motivated by their size, mechani- 
cal flexibility and chemical tunability. This vision is near- 
ing reality with the capacity to self-assemble, functional- 
ize, and reproducibly measure the current-voltage (I-V) 
characteristics of small groups of molecules. Molecular 
I-Vs have revealed a wide range of conductance proper- 
ties, from metallic conduction in carbon nanotubes |l]] 
and quantum point contacts J2j, to semiconducting be- 
havior in DNA H and conjugated aromatic thiols Q, and 
insulating behavior in alkyl thiol chains Q. Interesting 
device characteristics such as rectification ||, switching 
|tJ and transistor action || have also been reported. 

The classic paradigm for asymmetry in molecular I- 
V measurements is the Aviram-Ratner diode, consisting 
of a semi-insulating molecular species bridging an elec- 
tron donor-acceptor pair Q . A positive bias on the con- 
tact at the donor end brings the energy levels on the 
donor and acceptor sites into resonance, while the oppo- 
site bias moves the system away from resonance, leading 
to a strongly asymmetric I-V characteristic [|[@. Spa- 
tial asymmetry in the molecule is essential in generating 
the I-V asymmetry, causing the energy levels, the elec- 
trostatic potential and electron wavefunctions to be quite 
different for positive and negative voltages Jll| . 

In this paper, we address an I-V asymmetry exper- 
imentally observed for spatially symmetric molecules 
|L2| Q that is qualitatively different from and weaker 
than the rectification in spatially asymmetric molecules. 
The molecular I-V curves seen for these systems start off 
being symmetric, but pick up a weak, reversible asymme- 
try as the contacts are manipulated (Fig. 1). In contrast 
to spatially asymmetric molecules, conduction in these 
molecules at opposite voltages occur through essentially 
the same molecular levels with very similar wave func- 
tions. We show that the origin of the observed asymme- 
try is nontrivial, involving self-consistent shifts in the en- 
ergy levels due to charging effects. Asymmetry in charg- 
ing arises due to unequal coupling with the contacts, and 



can be present in conduction measurements performed 
with a break junctio n [fL2| , an STM tip or an evap- 
orated gold contact |lCf|. Remarkably, the sense of the 
asymmetry depends crucially on whether the conduct- 
ing level is HOMO or LUMO. We establish that for a 
spatially symmetric molecule the current is lower for pos- 
itive bias on the stronger contact if conduction is through 
a HOMO level, and higher if conduction is through a 
LUMO level. 
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FIG. 1. Sets of I-V characteristics (vertically shifted) ob- 
served by Reichert et al. |l2j for a break junction contacting 
a spatially symmetric molecule. To start with the I-V char- 
acteristics are symmetric (center). On weakening one of the 
contacts by pulling, the I-V curves become asymmetric. Two 
opposite sets of asymmetric I-Vs are obtained with the same 
molecule simply by pulling on the contacts (top and bottom) . 

Although ab-initio calculations exist for molecular con- 
duction p5|-p^], there is uncertainty about the position 
of the Fermi level even in idealized metal-molecular het- 
erostructures. This uncertainty translates into a corre- 
sponding ambiguity in the conductance gap in molecular 
I-Vs, as well as the nature of the conducting orbitals. For 
phenyl dithiol (PDT) and its derivative molecules cou- 
pled to gold (111) contacts, the Fermi level has variously 
been suggested to be mid-gap 18 1, closer to the HOMO 
level 0,|19],|O), or closer to the LUMO |||l|- Although 



1 



all these models predict similar symmetric I-Vs, the sense 
of the I-V asymmetry for asymmetric contacts pins down 
the nature of the conducting orbital rather restrictively; 
for instance, the I-V asymmetry measured for PDT with 
an STM tip H] indicates HOMO-based conduction, as 
we show later. 
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FIG. 2. Schematic description of resonant current conduc- 
tion through discrete molecular levels broadened by interac- 
tion with the contacts. Onset of current corresponds to cross- 
ing a molecular level by a contact chemical potential under 
bias. For a molecule strongly coupled to a substrate, and 
with a large gap DOS (a), the chemical potential of the other 
(weaker) contact crosses the HOMO and LUMO levels at op- 
posite bias. In a typical molecular conductor (b) on the other 
hand, the same level is crossed both ways. 

Non self- consistent description of molecular conduc- 
tion. An isolated molecule has discrete molecular lev- 
els, that are broadened upon coupling with contacts into 
a continuous density of states (DOS) due to hybridiza- 
tion with the metal wave functions. For molecules like 
PDT the broadening is weak enough that the DOS shows 
a clear HOMO-LUMO gap (HLG). The difference in the 
work functions between the metal and the molecule leads 
to charge transfer between them, and the molecular levels 
adjust due to Coulomb charging until the heterostructure 
is in equilibrium. Under bias the contact chemical poten- 
tials split by the applied voltage (/ia — Mi = eV) . As long 



as both potentials lie in the HLG, all HOMO levels are 
occupied and LUMO levels unoccupied, and there is es- 
sentially zero current. Once the bias is large enough that 
a contact chemical potential crosses a molecular level, 
that level is filled by one contact and emptied by the 
other, and starts conducting current. Thus the specific 
levels contributing to conduction for either bias direction 
depend on the way /xi and fi2 are disposed relative to the 
molecular levels under bias. 

One can get an asymmetric I-V with a symmetric 
molecule provided the latter is coupled strongly enough 
to the substrate to be in equilibrium with it, resulting 
in fi\ = Ep. The potential of the other contact (say, an 
STM tip) moves freely under bias, \i2 — Ep + eV (Fig. 
2a). If furthermore there is a large density of metal- 
induced gap states (MIGS) in the gap, then the electro- 
static potential and the molecular levels stay pinned to 
/ii and don't vary with V. For positive and negative 
substrate bias, /i2 crosses the LUMO and HOMO levels 
respectively, yielding a strongly asymmetric I-V, since 
the LUMO and HOMO transmissions are quite different 
in general. The conductance gap is then given by the 
HLG, a property intrinsic to the molecule itself. 

In practice, however, spatially symmetric molecular 
species typically exhibit symmetric I-Vs 0, with a con- 
ductance gap that varies from experiment to experiment 
jf].^2|. The MIGS are usually not appreciable caus- 
ing the electrostatic and chemical potentials to separate. 
The average electrostatic potential across the molecule 
is V/2, which shifts the levels along with the applied 
bias. Relative to the molecular levels therefore, the con- 
tact potentials move in opposite directions symmetri- 
cally, jLti,2 = Ef =p eV/2 (Fig. 2b). The same molecular 
level (HOMO in the figure) is crossed for either bias by 
the contact chemical potentials, leading to a symmetric I- 
V, with a conductance gap given by 4(_E/ — -Ehomo) (H) • 
The gap is no longer the HLG, but depends on both the 
molecular chemistry and the contact microstructure. 

The above non self-consistent description of molecular 
conduction can be put on a quantitative footing by us- 
ing an appropriate molecular Fock matrix F (ab-initio or 
semi-empirical), followed by a non-equilibrium Green's 
function (NEGF) formulation of transport For 
a given Fock matrix F, overlap matrix 5, and con- 
tact self-energies Xi : 2 with corresponding broadenings 

Ti.2 = i ^Si.2 — 2 ^ j the energy levels are given by the 
poles of the nonequilibrium Green's function G, while 
their occupancies are obtained from the corresponding 
density matrix p and the contact Fermi functions /j j; 

G{E) = {ES - F - Si - E 2 ) _1 

/OO 
dE (/iGTrGt + f 2 GT 2 G^) 
-OO 



hAE) - [1 + exp ((E - fn^/keTT 1 . 
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The number of electrons N and the steady-state current 
I are then given by: 

N = 2 (for spin) x trace(p S) 

Op f 00 

I = -^J dEtTa,ce(T 1 GT 2 G r )[f 1 (E)- f 2 (E)}. (2) 

The important point to notice is that since the same 
molecular level is being crossed both ways, the current 
is symmetric in the coupling constants Ti.2 and one is 
stuck with a symmetric I-V for a symmetric molecule, no 
matter how asymmetric the contacts are. 
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FIG. 3. Origin of asymmetry in charging: one side (say 
a substrate) is strongly contacted while the other is weak. 
Although the same level is crossed by the contact poten- 
tials for opposite bias, for positive substrate voltage (b) the 
HOMO level is emptied by the stronger contact, which posi- 
tively charges the molecule and shifts the energy levels down. 
Such a shift, not present for negative substrate bias (a), post- 
pones the onset of conduction and effectively stretches out the 
voltage axis in the direction of positive substrate bias. 

The asymmetric I-V cannot simply be explained by 
suggesting that the average electrostatic potential across 
the molecule is nV (^3|, n decreasing from 0.5 as a con- 
tact is drawn out. Such a theory yields an asymmetry 
only in the onset voltages for conduction, rather than a 
dragged out I-V as in Fig. 1 (the latter actually corre- 
sponds to different conductance peak values in the two 
bias directions). One will then need to invoke a voltage- 
dependent r) to explain Fig.l, which formally amounts to 
our self-consistent charging description below. 

Self- consistent description (Charging induced asymme- 
try) The above description of molecular conduction does 
not include the effects of charging and electronic Stark 
shift back on the levels. For a small molecule with a large 
(~ 1 eV) capacitative charging energy, charging effects 
need to be taken into account self-consistently. In the 
following, we show that such effects can lead to contact- 
induced asymmetries of the sort observed in Fig. 1. 

The effect of differential charging on the I-V asymme- 
try is explained schematically in Fig. 3 for a single level. 



Let us suppose that the contacts are strongly asymmet- 
ric (ri 3> T2). For negative bias on the strong contact 
(a), the latter is trying to fill the nearest (HOMO) level, 
while the weaker contact is trying to empty it, with the 
net result that the HOMO level stays filled, with cur- 
rent onset set by the voltage where \iv first crosses the 
(neutral) molecular level. For positive bias (b), however, 
the HOMO level is emptied out, which charges up the 
molecule positively. This adds a self-consistent charg- 
ing energy that lowers all the energy levels, postponing 
thereby the point where the HOMO is crossed by In 
effect, this stretches out the voltage axis, leading thus to 
a smaller conductance for positive bias on the stronger 
contact, as mentioned earlier. For LUMO-based conduc- 
tion the argument is reversed, since filling the LUMO 
level charges up the molecule negatively. 
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FIG. 4. Extended Hiickel-based I-V characteristics for 
Phenyl Dithiol (PDT), including self-consistent charging ef- 
fects. The central graph (symmetric) is for symmetric weak 
couplings (rs decreased symmetrically by a factor of 10 rel- 
ative to their chemisorbed values) of the PDT molecule with 
Au(lll) contacts. The lower graph corresponds to decreasing 
T2 five times relative to Ti, which stretches out the voltage 
axis in the direction of the stronger contact, i.e., along positive 
bias on contact 1. Reversing the couplings yields the upper 
graph. This graph can also be obtained without reversing the 
couplings, by simply raising Ef so that the conduction occurs 
through the LUMO level instead. 

We illustrate self-consistent charging effects in a PDT 
molecule modeled by a extended Hiickel |2j| description 
of the Fock matrix F, coupled to the NEGF equations 
described earlier. The self-energies are calculated for 
Au(lll) contacts using a recursive technique ^J§. Charg- 
ing effects are incorporated in F by adding a term Use 
describing the deviation from the equilibrium electron 
count Nq within the local density approximation 25 



U sc = U (N- N ) , (3) 
and solving for F self-consistently with equations [l] and ||. 
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The constant charging energy Uq (~ 0.7eV per electron) 
corresponds to a flat electrostatic potential profile in the 
molecule pjjj] , and serves to shift all the energy levels 
rigidly by a constant. 

Fig. 4 shows the I-V characteristic for PDT with weak 
symmetric contacts (central graph), with contact self- 
energies 0.1 times their chemisorbed values, assuming 
HOMO conduction (Ef = 10 eV). Weakening the cou- 
pling with the right contact (f^ = 0.2 Ti) stretches the 
voltage axis along contact 1 (lower figure) . The I-V asym- 
metry can be switched (upper figure) either by switching 
the contact couplings, or by raising Ef enough to gener- 
ate LUMO-based conduction. 

The calculated trend in the I-V characteristics 
agrees qualitatively with the experiments on symmetric 
molecules with asymmetric contacts [fL2|~[l4"|l . In break- 
junction experiments in general, the molecules are likely 
to adsorb separately onto the two contacts symmetrically 
p6j, leading to symmetric I-Vs even for asymmetric 
contacts. In carefully controlled break-junction experi- 
ments, however, a single molecule is expected to couple 
simultaneously with both contacts. This allows one to 
induce asymmetry easily by manipulating the contacts, 
changing the I-V from symmetric to asymmetric and 
back |l2j , as well as flipping the sense of the asymmetry 
by interchanging contact couplings. Weak I-V asymme- 
tries are also seen in STM measurements of symmetric 
molecules such as PDT [|l4j], where the conductance is 
lower for positive substrate bias. Since the coupling of 
the molecule with the STM contact is evidently weaker 
than with the substrate, the above asymmetry implies 
that conduction in PDT is HOMO-based. 

In summary, we have established that even for spatially 
symmetric molecules, contact asymmetry can induce an 
asymmetric I-V through differential charging, the sense 
of the asymmetry depending on whether conduction is 
through a HOMO or a LUMO level. 
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